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ABSTRACT 


The extinction paradox is a contradiction between results of geometrical optics which predict that 
at high frequencies the scattering cross section of an object should equal its geometrical cross 
section and rigorous scattering theory which shows that at high frequencies the scattering cross 
section approaches twice the geometrical cross section of the object, Confusion about the reason 
for this paradox persists today even though the nature of the paradox was correctly identified 
many years ago by Brillouin. In this paper the resolution of the paradox will be restated and 
illustrated with an example, and then the implications to the interpretation of scattering cross 


sections will be identified, 


Introduction 


The scattering cross section, a s , for a particle is defined to be the ratio of the total power 
scattered by the object to the power per unit area incident on the object (e.g, Bom and Wolf, 
1959; Bohren and Huffman, 1983; Ishimaru, 1978; van de Hulst, 1981): 

total power scattered 

= ( 1 ) 

incident power per unit area 

To calculate the scattered power it is standard procedure in electromagnetic theory to resolve the 
total electric and magnetic fields into “incident” and “scattered” components as follows: 

^ “ ^inc * ^scat (-) 

^ ” ^inc ^ ^scat 

Then the time average “scattered” power density can be obtained from the Poynting vector for 
the scattered fields: 


scat 


(1/2) Re [E scat 


(3) 


and the scattering cross section becomes: 


ff 


scat 


ds 


Pine 


(4) 


A. 


where i is a unit vector in the direction of propagation of the incident wave, Pj nc = (1/2) Re 
[Ej nc x H£, c ] and the integration of Equation 4 is over a closed surface surrounding the particle. 


In the case of particles which are large compared to wavelength (i.e, in the limit k-*») a simple 
intuitive solution for a s can be obtained from geometrical optics. To do so consider a plane wave 
with wave number k = oj ,//ae"incident on an object of finite size as illustrated in Figure I. In the 
limit k-*», one expects the object to cast a shadow on a screen placed behind the object, and one 
expects this shadow to be approximately the orthogonal projection of the object onto the screen. 
Such a shadow is called the geometrical cross section of the object. If the object is completely 
opaque, then the shadow is a region in which there are no electromagnetic waves and therefore 
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no power. Under these circumstances, the object must have scattered (or possibly absorbed) an 
amount of power equivalent to that which in the absence of the particle would have been in the 
shadow. This scattered power can be calculated by multiplying the Poynting vector of the incident 
wave times the geometrical cross section of the object. Hence, if no power is absorbed by the object, 
then the scattering cross section should be equal to the geometrical cross section of the object. (If 
the object is absorbing, then the scattering cross section will be smaller than the geometrical cross 
section depending on the amount absorbed.) 

The preceding derivation of the scattering cross section is based on intuition obtained from 
geometrical optics. It is correct to the same extent that geometrical optics is correct and therefore 
ought to apply when k->«». However, the expression obtained for the scattering cross section is not 
consistent with solutions obtained using exact solutions for the scattered fields, and H scat . 

This was apparently first noticed in connection with the Mie solution for scattering from a perfectly 
conducting sphere (Brillouin, 1949). The geometrical cross section of a sphere of radius a is 7ra 2 but 
the total power scattered from the sphere divided by the incident power density approaches 27ta 2 in 
the limit k-*» (Born and Wolf, 1959; van de Hulst, 1981). The same is true of other “canonical” 
shapes for which exact solutions to the wave equation can be found. For example, the geometrical 
cross section per unit length of a perfectly conducting infinite circular cylinder of radius a is 2a, but 
the scattering cross section per unit length obtained from theory approaches 4a in the limit k-*» 
(Bowman, Senior and Uslenghi, 1969). The scattering cross section also approaches twice the geo- 
metrical cross section for scattering from cylinders of elliptical cross section, independent of polari- 
zation, and this is also true for acoustic scattering from spheres (Bowman. Senior, and Uslenghi, 1969) 

Hence, theory predicts that the scattering cross section approaches twice the geometrical cross 
section as k-*» whereas geometrical optics arguments lead to the conclusion that the scattering cross 
section is equal to the geometrical cross section. This contradiction is not dependent upon assump- 
tions made about absorption because the results cited above are for lossless (non-absorbing) objects. 


Nor is it the result of approximations inherent in geometrical optics (e.g. which neglects dif- 
fraction). In fact, both arguments described above are correct in the limit k-*». The contra- 
diction is of more subtle origin and has been called the “extinction” paradox (e.g. Ishimaru, 

1978; Bohren and Huffman, 1983; van de Hulst, 1981.) 

The extinction paradox is a consequence of the definition (Equation 2) adopted for the 
scattered fields. This was originally pointed out by Brillouin (1949) and is particularly evident 
in the limit k-*». In this case, E scat and R scat correctly give the fields reflected from the object 
but as defined in Equation 2 they also include a term which combines with the incident wave to 
produce the shadow behind the object. For example, consider a plane mirror in the presence of 
radiation from a high frequency source as illustrated in Figure 2. In the limit k-*» one expects 
to find a reflected wave and also a shadow behind the mirror. Thus, E scat as defined in Equa- 
tion 2 must consist of one term which combines with the incident field to produce the shadow 
and must also contain a term which gives the wave reflected by the mirror. In the limiting case 
(k-*») there is no power in the region behind the mirror (e.g. a completely dark shadow) and the 
mirror simply redirects the incident power into the reflected wav-:. The power in the reflected 
wave is therefore just the power intercepted by the mirror from the incident wave, and in this 
case the ratio of the reflected power divided by the incident power per unit area is just the geo- 
metrical cross section of the mirror, This is the geometrical optics argument for deriving the 
scattering cross section. However, when one computes the total scattered power using Equation 
3, one is including in E^ and H^ a contribution which accounts for the reflected wave and also 
a contribution from the terms needed to create the shadow. Clearly, in the limiting case (k-*») 
with a completely opaque, lossless mirror, these two contributions are equal and each is equal to 
the power intercepted from the incident wave by the geometrical cross section of the object. Thus, 
in the high frequency limit one expects the scattering cross section obtained from Equation 4 to 
yield twice the geometrical cross section, as in fact is found to be true when the calculation is 
performed for simple shapes using the “exact” solutions for E scat and H scat , 
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Example 


The two component nature of the scattered fields will be illustrated in this section with an 
example of scattering from dielectric disks. Imagine a plane wave incident on a dielectric disk as 
illustrated in Figure 3, The unit vector, I s , in the direction of propagation of the incident wave 
lies in the y-z plane and the unit vectors 'll and “v, indicating the direction of the electric field for 
horizontally and vertically polarized waves, are perpendicular or parallel to the y-z plane, respec- 
tively. The unit vector o indicates the direction of propagation of the scattered wave. Calculations 
of the scattered electric field, E scat , due to the incident plane wave have been made using the 
physical optics approximation of Le Vine, et al., (1983). In this approach, a solution for the 
scattered fields is first written in terms of the sources (currents) induced in the disk by the inci- 
dent wave. These are unknown, but for linear media can be expressed in terms of the dielectric 
properties of the medium (complex dielectric constant) and the " .Ids inside the disk, The approx- 
imation is to use the fields inside an identically oriented, ■nfimi, 1 ' slab of the same thickness and 
dielectric properties to approximate the actual fields inside the disk. One expects the results to 
be reasonable for non-grazing incidence for disks whose minimum characteristic dimension (e.g. 
diameter of circular disks) is very large compared to wavelength. It has been shown (Le Vine, 

1984) that this approximation is consistent with results obtained using other techniques and 
agrees with experiments at normal incidence. 


To perform the required calculations it is convenient to express the scattered fields in terms 

r ,AAx 


of a dyadic scattering amplitude, f (o, l), as follows: 

Escat * P • f(o 


A =,A A eik r 


(5) 


A . 


where p is a unit vector in the direction of polarization of the incident wave and E 0 is the ampli- 
tude of the incident wav* , Thus, the component of the scattered wave with polarization in the 
q-direction due to an incident wave with polarization p is: 


p A A = A A A e jkr 

E scat * 9 = P * f (o, l) • q — fc 0 


( 6 ) 
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fit 


Calculations have been made of'v • f (o,1) * v using the physical optics theory of Le Vine, et al, 
(1983; Equations 22-23), Examples are shown in Figure 4 for a circular disk of radius 10 cm, 
thickness 0.5 cm, frequency 9 GHz (ka = 6rr) and relative dielectric constant of e r * 25 + jl 1 
which is representative of leaves at microwave frequencies, Polar plots of | v • f (o,^ • v'i as a 
function of o are presented for fixed 'll Each plot was obtained by calculating (v • Ffo,^ • v| as 
o was rotated through 360° in the z-y plane, The z-axis (o=i) is at the top of each figure and 
the circles are contours of constant amplitude, the outer most circle representing unit amplitude 
and the inner most circle being 0,2, Four plots are shown in Figure 4, They differ only in the 
angle of incidence of the plane wave (which is indicated by the arrow), In the example labelled 
0 = 0°, the plane wave is normally incident on the disk (f=-n = -z) and in successive examples f 
is obliquely incident at increasingly large incidence angles, 0. Notice that in each example the 
magnitude of the scattering amplitude has two peaks: one which is always behind the disk in the 
direction of propagation of the incident wave, and another which is in the direction of specular 
reflection (angle of incidence equals angle of reflection). The peak in the forward direction is 
due to the component of the scattered field which combines with the incident wave to form the 
“shadow” behind the disk, and the peak in the direction of specular reflection represents energy 
whose direction of propagation has been changed by the disk (i,e. electromagnetic fields reflected 
from a mirror), The two-peaked nature of E scat evident in Figure 4 is not peculiar to disks nor to 
the method used here to obtain a solution; rather it is characteristic of scattering at high fre- 
quencies (e,g, Papayiannakis and Kriezis, 1983, and Kozaki, 1983). 

The two peaks in the scattered field both contribute to the “scattered” power and therefore 
to the scattering cross section of the disk. In particular, the power density with polarization v 1 
scattered from the disk in the o direction by an incident wave with polarization 'v is : 

F W = Vi Re (E scat x H* cat • o) (7) 

= Vt y/ejjl Eq |v • F (o.r) • v| 2 
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where Vi v/e/jx Eq is the power density in the incident wave, 

ah directions divided by the incident power density is 

Total power scattered with polarization v 
Incident power density with polarization y 

It is clear from the example in Figure 4 that the integral in Equation 8 includes contributions from 

both the fields specularly reflected from the disk and also from the fields needed to create the 

shadow. 

IMPLICATIONS 

As indicated in the previous paragraphs, the extinction paradox is the result of a contradiction 
between physical intuition which leads one to understand the term “scattered’’ to mean “reflected”, 
and the conventional definition of scattered fields which includes both the reflected fields and a 
term which combines with the incident wave to create a shadow behind the object. Problems 
occur when one uses the conventional definition but interprets the results using the intuitive 
notion. Two examples will be considered below. 

A. The extinction cross section 

The time average power absorbed by a particle, P a , is equal to the net flux of the time average 
Poynting vector into a dosed surface surrounding the partide. That is: 

-P a = H Re ff [Einc + E’scat 1 * [H* c + H* at ] • ds (9a) 

= Ps -Pt 

where 

p s = M Re ff %cat X f! scat * d5 W 

- P t = V 2 Re ff [E^c x H* at + ^ cat x Hf nc ) • dS (9c) 

Dividing both sides by the incident power density (i.e. P^ = Vi s]7\\i Eq ) and using the definitions 
°a = Pa/Pinc» a s = p s / p inc and °T = p t /Pinc> one may write Equation 9a in the form: 

o T = a s + o a ( 10 ) 


Thus, the total power scattered in 
JJ l v • F (0,1^) ♦ v| 2 d£2 (8) 
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where op is called the “total” or “extinction” cross section and o a and a s are called the “absorp- 
tion and “scattering” cross sections, respectively. By substituting Equation 5 into Equation 9c one 

A 

can express the total cross section associated with an incident wave with polarization p in terms of 
the scattering amplitude as follows (e.g. Bom and Wolf, 1959; Karam and Fung, 1982): 

oj - (2?r/k) Im p • F(i^ • p (1 1) 

This result is called the “optical” or “forward” scatter theorem. 

The conventional interpretations of ay is as a measure of the power scattered and absorbed by 
the object (e.g, Bohren and Huffman, 1983; van de Hulst, 1981; Ishimaru, 1978); specifically, that 
it is an equivalent area which intercepts from the incident wave the same power that is absorbed 
and scattered by the particle. This interpretation would be correct if o s were a measure of the 
power “reflected” from the object, However, as was pointed out above, cr s also includes a contri- 
bution from the shadow producing fields. Although these terms may not be distinct for finite k. 
one can write for conceptual purposes: 

CT s = °rs + a ss 

where a rs accounts for the reflected power and o ss is the contribution due to the shadow producing 
part of E scat . Using this notation one can write: 

ay = (cr ts + % ) + a a 

= (o rs + o a ) + o ss (12) 

= °t + <%s 

where a t = o rs + a r Clearly, o t is the cross section which measures the t yower reflected and ab- 
sorbed by the particle, and ay measures the power reflected and absorbed plus an additional 
term which is present because, by nature of their definition, ^ cat and Hs cat include a term which 
combines witli the incident fields to produce a shadow behind the object. When there is no ab- 
sorption (o a = 0) and the object is completely opaque, conservation of energy leads one to expect 
that o rs = o ss : and when k-*<» geometrical optics suggests that a rs = o g where a g is the geometrical 

cross section of the object. Hence for lossless, opaque objects one expects: 

*7 


Jim o t » 2o b 


(13a) 


lim o t « o g 
k-**> 


Cl 3b) 


Thus, the definitions are consistent with both intuition and theory. Paradoxical results occur 
when one tries to interpret oj as a measure of the total power "reflected” and absorbed without 
consideration of the fact that it includes an additional term which is present as a consequence of 
the definition adopted for ^scat ■ In fact, Equations 10 and 11 are correct, but the label "total” 
or "extinction” cross section given to oj is misleading, 


B. Radar Cross Section 

It is conventional to define the radii*' cross section of an object as follows (Bowman, Senior 
and Uslenghi, 1969; Ishimaru, 1978; Ruck, et al., 1970): 


o - lim 4?rR 2 

R->oo 


E 


scat 


E* 


scat 


(14) 


I ^inc * Eg lc I 

In the case of an incident wave with polarization p and a scattered wave with polarization q one 
obtains (e.g, using Equation 5): 


a pq = 4tr| p • F(o,'i} • q | 2 (15) 

The physical interpretation of Equations 14 and 15 is that o p q is 4 rr times the power per unit solid 
angle scattered in the o direction normalized by the incident power density (e.g. see Equation 7). 
This is correct if E^ represents the scattered (i.e, reflected) power, But as pointed out above, 
E sca t a ^ so includes the shadow producing term. 


To illustrate the problems that can be encountered, again consider scattering from a disk at 

high frequencies. The scattering cross section o v w behaves as a function of o as indicated in 

Figure 4. (See Le Vine, 1984 for examples in a more conventional format). Because the frequency 

is high (large ka) c vv has two distinct peaks which are manifestations of the reflection and shadow 

producing terms in E sca t. At normal incidence it can be seen from Figure 4 that the predominate 
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contribution to the backscattcr radar cross section (o is the reflected term in 'E scal . In this 
case one expects a reasonable agreement between the definition (Equation 14) and observation 
(i.e. a measurement) because the peak is, in fact, due to '‘reflected" power, However, in the case 
of forward scatter, the radar cross section as given by Equation 14 will also have a peak, However, 
this peak is due to the shadow producing term in E scat and does not represent “reflected” power. 
In the laboratory a radar receiver directly behind the disk will not receive the electric field, E scat ; 
rather it will be in the shadow of the disk where it will receive only a small signal. It is necessary 
to use E « E^ + E^ in this region to correctly describe the electric field at the receiver. 

Summary 

As has been pointed out above, the electric field, E scat , as conventionally defined (E - Ej nc 
+ E^) includes both waves “reflected” from the object and also a term which when combined with 
Ej nc produces the shadow behind the object, Thus, E scat is a measure of more than just the 
electric fields removed from the incident wave by the object. Failure to recognize this leads to 
inconsistencies such as the extinction paradox and to misleading interpretations such as that the 
total scattering cross section, <7j * a a + o s> is proportional to the power absorbed and scattered 
by the particle. On the other hand, as long as one uses both the incident and scattered fields (i.e. 
Ejnc + E sca t ) to represent tire total electric field, then the cross sections which occur in an analysis 
will be used correctly. For example, consider the propagation of a plane wave with polarization 
p through a layer of random scatterers. Using a Born approximation, one concludes (e,g. Bohren 
and Huffman, 1983; Oguchi, 1983) that the effective propagation constant in the layer is: 

k = kg + (2jr/ko) N (p • Fdl'i) • p) 

where N is the number of particles per unit volume, Thus, the effective attenuation coefficient, 
a = Im(k), for propagation through the layer is; 

a = N Im t(27r/k 0 ) p* f • p] 
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Tills result, which is derived using E ■ E^ + T^ at is correct to within the limits set by the Born 
approximation, What is incorrect (perhaps the more appropriate word is “misleading”) is to 
interpret this result to mean that because aj a o t + a y the attenuation is determined by N times 
the power “scattered” and absorbed in the medium. The attenuation coefficient is correctly given 
by N aj, but the scattering cross section, a v does not correctly measure the power deflected by 
the particles from the incident plane wave. In fact, in the high frequency limit it is incorrect by 
a factor of two (a s r, 2 a g ). 
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